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Understanding  the  lithium-silicon  alloying  behavior  is  essential  for  achieving  maximum  charge  capacity 
in  the  negative  electrodes  of  lithium-ion  batteries.  Our  atomistic  simulations  show  that  in  amorphous 
silicon  with  a  disordered  network  structure,  inserted  lithium  atoms  can  find  equilibrium  positions  in  the 
interstices  of  big  rings.  Alternatively,  lithium  is  incorporated  into  the  network  by  the  destruction  and 
reformation  of  smaller  rings.  These  atomic-level  mechanisms  are  characterized  by  using  the  network 
topology  measure  of  ring  statistics,  which  are  correlated  to  the  lithiation  responses  of  silicon  electrodes. 
The  results  reveal  the  influence  of  lithium  concentrations  on  the  electro-chemical-mechanical  behavior 
of  silicon.  Implications  on  the  reversibility  and  dynamics  of  the  lithiation  process  are  discussed. 

©  2010  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Rechargeable  lithium-ion  batteries  are  currently  the  preferred 
form  of  electrical  energy  storage  devices  [1,2].  They  can  deliver  high 
energy  density  (~210Whkg_1;  650  Whl-1)  with  long  life  cycle. 
Their  recent  and  potential  applications  include  portable  electronic 
devices,  cordless  power  tools,  and  electric  vehicles.  Silicon  is  being 
considered  as  an  anode  material  for  Li-ion  batteries  because  it 
has  the  highest  known  theoretical  charge  capacity  (4200  mAh  g-1 ) 
[3-8].  However,  Si  anodes  often  suffer  from  pulverization  and 
capacity  fading.  This  is  caused  by  the  large  volume  changes  of 
Si  (~300%)  upon  Li  insertion/extraction  close  to  its  theoretical 
charge/discharge  limit.  This  issue  is  particularly  significant  when 
the  Si  electrodes  are  subject  to  high-rate  charging  or  discharging, 
which  is  highly  desirable  but  often  induces  a  non-uniform  distri¬ 
bution  of  Li  in  Si.  The  resulting  large  incompatible  deformation 
between  areas  of  different  Li  contents  tends  to  initiate  fracture, 
leading  to  the  electro-chemical-mechanical  failure  of  Si  electrodes 
[9-11].  Are  there  new  ways  to  optimize  the  charge  capacity  of 
Si  anodes  while  preventing  failures?  To  address  this  question, 
understanding  the  Li-Si  alloying  behavior  is  essential,  requiring 
fundamental  research  that  goes  beyond  empirical  battery  cell  stud¬ 
ies. 
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Metastable  amorphous  glass  is  formed  when  the  crystallization 
of  intermediate  compounds  is  frustrated  [4].  When  reacted  with 
Li,  Si  crystals  yield  solid-state  amorphous  phases  at  low  tempera¬ 
ture,  but  crystalline  phases  at  high  temperature  (415  °C)  [12].  The 
formation  of  amorphous  LixSi  alloys,  where  x  denotes  the  ratio  of 
Li-Si,  has  been  experimentally  observed  at  room-temperature  in 
the  electrochemical  reactions  of  Li  with  Si  thin  films  [4],  nanoparti¬ 
cles  [5]  and  nanowires  [6].  In  principle,  amorphous  LixSi  avoids  the 
formation  of  coherent  boundaries  between  coexisting  crystalline 
phases  of  different  Li  concentrations.  This  can  mitigate  the  buildup 
of  incompatible  strains  (and  thus  internal  stresses)  in  Si  electrodes 
so  as  to  improve  the  cycling  performance  of  batteries.  In  addition, 
a  variety  of  Si  nanostructure  morphologies  have  been  explored  to 
enable  Si  to  perform  well  as  anodes  [5-7],  as  it  is  supposed  that  the 
nanometer  size  scale  can  facilitate  strain  relaxation,  enhance  flaw 
tolerance  [13],  shorten  diffusion  path,  and  increase  surface  area  of 
Si  to  better  react  with  Li. 

As  a  fundamental  step  towards  understanding  the  Li-Si  alloy¬ 
ing  behavior,  here  we  report  an  atomic-level  study  of  Li  insertion 
in  amorphous  Si  (a-Si).  The  structures  and  energetics  of  Li-Si 
alloys  are  computed  quantum-mechanically  by  using  the  molec¬ 
ular  orbital  theory  [14].  Detailed  analyses  of  the  ring  geometry 
and  statistics  are  related  to  various  electro-chemical-mechanical 
characteristics  of  Li  reactions  with  a-Si.  They  are  found  to  depend 
sensitively  on  the  local  Li  concentration.  The  results  provide  an 
atomistic  basis  for  further  study  of  the  dynamic  alloying  processes 
and  electrochemically  driven  solid  state  amorphization  between  Li 
and  Si. 
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Fig.  1.  Comparison  of  the  amorphous  and  diamond  cubic  crystals  of  Si  based  on  radial  and  angular  distribution  functions  (RDFs  and  ADFs). 


2.  Methods 

The  a- Si  structures  are  generated  from  diamond-cubic  Si  crys¬ 
tal  by  using  an  activation-relaxation  technique  (ART)  [15].  As  a 
computationally  efficient  approach,  ART  focuses  on  simulating  the 
activated  events  of  escaping  the  basins  of  attraction  in  the  energy 
landscape,  i.e.,  they  jump  over  the  barriers  that  separate  different 
local  energy  minima.  ART  can  be  used  to  avoid  the  bias  on  struc¬ 
tural  evolution  from  the  fast  cooling  rate  in  the  melt-and-quench 
molecular  dynamics  method  [  1 6].  To  speed  up  calculations,  the  a-Si 
structures  are  first  created  by  using  the  empirical  Stillinger-Weber 
[17]  interatomic  potential  and  then  further  relaxed  by  the  Molec¬ 
ular  Orbital  PACkage  (MOPAC  2009)  with  the  AMI  method. 

We  adopt  the  protocols  developed  by  Chevrier  and  Dahn  [18,19] 
for  modeling  the  lithium  insertion  and  extraction.  The  Li  atoms  are 
randomly  inserted  into  a-Si  one  at  a  time  until  x  =  4.  After  the  inser¬ 
tion  of  each  Li  atom,  the  system  is  fully  relaxed  under  periodic 
boundary  conditions;  this  involves  the  relaxation  of  the  internal 
atomic  coordinate,  supercell  size  and  shape.  The  lithium  extrac¬ 
tion  process  is  similarly  modeled.  Starting  from  pure  a-Si,  we  have 
repeated  the  insertion-extraction  cycle  several  times,  and  found 
that  the  electro-mechanical  responses  of  volume  expansion  and 
electrical  potential  will  reach  a  steady  state  after  about  three  cycles. 
To  determine  the  electrical  potential,  we  note  that  neglecting  the 
entropy  contributions  [18,19],  the  free  energy  of  the  system  can 
be  approximated  by  the  potential  energy  which  is  a  function  of 
atomic  coordinates  at  zero  temperature.  Then  the  electrical  poten¬ 
tial  of  the  LixSi  structure  versus  Li/Li+,  in  units  of  electron  volts,  is 
calculated  by  V  =  -( dEhixSi/dx  -£Li),  where  £LixSi  is  the  potential 
energy  of  the  LixSi  structure  divided  by  the  number  of  Si  atoms 
in  the  system  and  ELi  is  the  energy  of  a  single  Li  atom  in  amor¬ 
phous  Li.  To  obtain  the  average  responses  for  comparison  with 
experimental  measurements,  the  potential  and  volume  are  aver¬ 
aged  over  the  calculation  results  of  six  statistically  independent  a-Si 
structures  at  the  respective  steady  states.  The  averaging  method 
that  we  used  is  the  Savitzky-Golay  smoothing  filter  with  the  main 
advantage  of  preserving  features  of  the  distribution  such  as  relative 
maxima,  minima  and  width,  which  are  usually  ‘flattened’  by  other 
techniques  (like  moving  averages,  for  example)  [20]. 

3.  Results  and  discussion 

The  starting  amorphous  structure  in  our  calculations  has  20  Si 
atoms.  The  supercell  volume  is  428  A3  corresponding  to  a  typical 
density  of  a-Si,  2212  kg  m-3  [21].  While  one  would  like  to  use  a 
larger  starting  system,  we  note  that  after  Li  insertions,  the  amor¬ 
phous  LixSi  alloy  at  its  charging  limit  involves  about  four  times 


more  atoms.  This  limits  the  size  of  the  starting  state  because  of 
the  computational  cost  associated  with  a  large  number  of  calcula¬ 
tions,  which  are  necessary  for  averaging  the  voltage-composition 
responses  and  for  analyzing  the  statistics  of  amorphous  states.  To 
validate  the  use  of  such  a  small  starting  system,  we  have  generated 
a-Si  with  various  sizes,  and  verified  that  it  well  represents  the  bulk 
behavior  in  terms  of  providing  similar  structural  characteristics  of 
the  short-  and  medium-range  order.  The  former  can  be  character¬ 
ized  by  radial  and  angular  distribution  functions  (RDFs  and  ADFs), 
and  the  latter  by  ring  statistics.  As  shown  in  Fig.  1 ,  RDF  has  a  narrow 
nearest-neighbor  peak,  while  the  peak  of  ADF  centers  around  109°, 
characteristic  of  tetrahedral  Si  bonding.  These  features  in  RDF  and 
ADF,  which  are  consistent  with  previous  modeling  results  [21,22], 
indicate  that  the  local  atomic  environments  of  a-Si  are  similar  to 
those  of  atoms  in  the  crystalline  phase,  as  typically  found  in  the 
solid-like  amorphous  covalent  structure  with  a  continuous  random 
network  [17]. 

The  incorporation  of  Li  into  a-Si  induces  two  kinds  of  struc¬ 
tural  changes  that  dominate  at  low  and  high  Li  concentrations, 
respectively:  insertion  in  the  interstitial  hole  of  big  rings  and  incor¬ 
poration  into  the  network  by  the  destruction  and  reformation  of 
small  rings.  We  find  that  when  the  Li  to  Si  ratio,  x,  is  less  than  about 
0.8,  the  added  Li  typically  resides  in  the  centers  of  the  big  rings 
after  structural  relaxation.  This  can  be  seen  in  Fig.  2(a)  and  (b), 
showing  an  8-membered  ring  of  Si  (highlighted  in  green)  before 
and  after  the  insertion  of  one  Li  (in  red),  respectively.  In  this  case, 
the  big  Si  ring  encloses  a  large  open  space  and  its  center  provides  an 
energetically  favorable  interstitial  site.  After  Li  insertion,  the  Si-Si 
bonds  in  the  ring  belt  do  not  break,  but  the  Li  atom  forms  new 
bonds  with  the  surrounding  Si  atoms  in  the  belt,  causing  a  small 
adjustment  of  the  ring  geometry.  Furthermore,  the  addition  of  this 
Li  divides  the  big  8-membered  ring  into  several  smaller  rings  (three 
3-membered,  one  4-membered  and  one  5-membered  rings).  Con¬ 
sequently,  the  fraction  of  big  rings  in  the  system  decreases  with 
increasing  x. 

Fig.  2(c)  and  (d)  shows  the  representative  structures  before  and 
after  the  addition  of  one  Li  atom  (in  red)  when  x>  0.8.  Because  of 
the  paucity  of  big  rings,  the  added  Li  has  to  typically  react  with 
surrounding  small  rings.  It  pushes  the  nearby  atoms  away  so  as  to 
create  a  larger  local  free  volume  for  its  incorporation  into  the  amor¬ 
phous  network.  In  conjunction  with  large  atomic  displacements  of 
surrounding  atoms,  the  local  ring  structure  is  reconfigured,  involv¬ 
ing  both  bond  breaking  and  formation.  For  example,  prior  to  the  Li 
addition,  the  length  of  a  homopolar  bond  between  Li  atoms  e  and  b 
in  Fig.  2(c)  is  2.69  A.  It  is  close  to  the  equilibrium  Li— Li  bond  length 
of  2.65  A  in  the  crystalline  phase  of  Lii5S4  from  density  functional 
theory  calculations  [23].  With  the  addition  of  a  Li  atom  (red)  in 
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Fig.  2.  Atomic  structures  of  lithiated  a-Si.  Large  and  small  spheres  represent  Si  and  Li,  respectively,  (a)  and  (b)  Before  and  after  the  interstitial  insertion  of  a  Li  atom  (red)  in 
an  8-membered  ring  of  Si  whenx  =  0.  (c)  and  (d)  Incorporation  of  Li  into  the  network  through  breaking  bonds  (in  red)  and  reformation  when  x  =  1.2.  (For  interpretation  of  the 
references  to  colour  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  the  article.) 


Fig.  2(d),  both  Li  atoms  e  and  b  form  new  homopolar  bonds  with 
the  inserted  Li,  while  the  bond  between  b  and  e  is  broken,  with  the 
separation  increased  to  4.05  A.  Evidently,  the  impact  of  Li  insertion 
on  the  local  ring  structure  involves  the  destruction  of  one  or  a  few 
small  rings  with  the  concomitant  reformation  of  several  similar 
sized  (3  or  4-membered)  small  rings. 

The  foregoing  atomic-level  structural  changes  associated  with 
Li  insertion  can  be  correlated  to  the  evolution  of  ring  statistics, 
which  measure  the  medium-range  order  of  amorphous  network 


structures.  We  use  Dijkstra’s  shortest  path  algorithm  [24],  imple¬ 
mented  in  MATLAB,  and  the  ring  counting  method  developed  by 
Franzblau  [25].  As  shown  in  Fig.  3(a),  the  ring  fraction  in  the  initial 
structure  of  a-Si  (x  =  0)  exhibits  a  single  peak  at  the  6-membered 
ring.  Similar  to  the  8-membered  ring  shown  in  Fig.  2,  this  kind 
of  big  6-membered  ring  provides  a  large  open  space  to  accom¬ 
modate  the  interstitially  inserted  Li.  Fig.  3(b)  indicates  that  with 
increasing  x,  a  peak  of  small  3-membered  rings  develops,  while  that 
of  6-membered  rings  diminishes.  This  change  is  understandable 


Ring  size  (x=0) 


Ring  size  (x=0.75) 


Fig.  3.  Histogram  of  the  distribution  of  ring  size  (i.e.,  the  number  of  members  in  the  ring)  in  lithiated  a-Si  with  different  ratios  of  Li  to  Si,  x.  The  fitting  curve  (in  red)  is  plotted 
to  guide  eyes.  (For  interpretation  of  the  references  to  colour  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  the  article.) 
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Fig.  4.  Comparison  of  the  theoretical  predictions  and  experimental  measurements 
[27,28]  of  the  electro-mechanical  response  of  lithiated  a-Si.  (a)  Ratio  of  the  current 
and  initial  volume  versus  Li  fraction  x;  (b)  voltage  versus  Li  fraction.  The  red  solid 
line  represents  the  statistical  average  using  the  Savitzky-Golay  smoothing  filter  [20] 
as  discussed  in  Section  2.  (For  interpretation  of  the  references  to  colour  in  this  figure 
legend,  the  reader  is  referred  to  the  web  version  of  the  article.) 

by  noting  the  earlier  observation  that  interstitially  added  Li  often 
divides  the  big  ring  into  several  smaller  rings.  Moreover,  we  show  in 
Fig.  3(c)  that  when  x  is  increased  to  around  0.8,  the  system  consists 
dominantly  of  3-membered  rings  with  a  single  peak.  Recall  that  in 
Fig.  2(d),  a  further  addition  of  Li  causes  the  reconstruction  of  local 
ring  structures  whenx  >  0.8.  This  can  be  correlated  to  Fig.  3(d)  show¬ 
ing  an  increase  of  the  single  peak  of  3-membered  rings.  Finally, 
we  note  that  a  small  percentage  of  large  rings  always  retain  even 
when  x  approaches  its  theoretical  limit.  This  arises  because  some 
large  rings  counted  in  the  statistics  correspond  to  the  case  of  small 
rings  running  along  the  belt  of  the  big  ring  [26].  These  small  rings 
occupy  the  open  space  enclosed  by  the  ring  belt.  As  a  result,  such 
large  rings  could  survive  because  there  is  little  room  to  facilitate  Li 
insertion. 

Fig.  4  shows  the  calculated  volume  and  voltage  versus 
Li  composition  curves,  which  are  in  accord  with  experimen¬ 
tal  measurements  of  lithiated  a-Si  [27,28].  The  slope  of  the 
volume-composition  curve  corresponds  to  the  volume  increase  per 
Li  insertion,  equivalent  to  the  so-called  partial  molar  volume  of  Li, 
Q.  Incidentally,  £2  is  a  key  parameter  in  the  coarse-grained  con¬ 
tinuum  model  of  diffusion-induced  stress  and  cracking  [9-11].  As 
seen  from  Fig.  4,  Q  is  small  and  non-linear  (the  numerical  fitting 
gives  £2  =  (0.57x2  -  0.04x+ 1.1 8)A3  per  Li)  when  x<0.8  (indicated 
by  the  green  dashed  line),  while  £2  is  large  and  nearly  a  constant 
(if2  ^  17.4  A3  atom-1  «1.1  x  10-5  m3  mol-1)  as  x>0.8.  This  trend 
can  be  correlated  to  the  foregoing  analyses  revealing  two  kinds  of 
structural  changes  that  dominate  at  low  and  high  Li  concentrations, 
respectively:  Li  interstitial  insertion  and  incorporation  into  the  net¬ 
work.  In  the  former  case,  the  interstitial  hole  can  well  accommodate 
the  inserted  Li,  so  that  the  induced  volume  expansion  is  relatively 
small.  This  is  in  contrast  to  the  latter  case  where  the  Li  incorporation 
occurs  through  network  reconstruction,  causing  a  larger  volume 
expansion.  Furthermore,  the  higher  voltage  at  small  x  is  related  to 
the  dominant  response  of  bond  formation  between  Li  and  Si  in  the 
ring  belt,  while  the  reduced  voltage  at  larger  x  is  attributed  to  the 


increased  energy  cost  associated  with  bond  breaking  and  stretching 
manifested  as  the  increased  volume  expansion. 

4.  Concluding  remarks 

In  summary,  the  present  study  reveals  the  atomic  structure, 
volume  and  voltage  responses  of  lithiated  a-Si.  The  results  rein¬ 
force  the  notion  that  the  Li-Si  alloying  behavior  and  associated 
properties,  including  the  network  topology  characteristics  of  ring 
statistics,  are  sensitive  to  the  local  Li  concentration  [16].  Our  work 
provides  an  atomistic  rational  of  the  relatively  small  volume  expan¬ 
sion  when  the  lithium  fraction  is  less  than  about  0.8.  This  effect 
arises  because  of  the  limited  dilatational  effect  associated  with 
the  insertion  of  lithium  atoms  that  mostly  reside  in  the  ener¬ 
getically  favorable  sites  of  interstitial  holes  in  the  open  network 
structure  of  Si.  Since  such  small  volume  expansion  has  an  intrinsic 
atomistic  structural  basis,  it  may  have  implications  for  reduc¬ 
ing  the  stress  buildup  and  thereby  improving  the  reversibility  of 
charging  and  discharging  processes  by  controlling  the  extent  of 
lithiation. 

While  the  work  focuses  on  the  static  aspects  of  lithiation,  it  also 
has  implications  on  the  dynamic  processes  of  Li-Si  mixing  and 
electrochemically  driven  solid  state  amorphization.  Consider,  as 
an  example,  Li  diffusion  in  a-Si  that  underlies  the  aforementioned 
processes.  It  is  generally  recognized  that  diffusion  in  amorphous 
materials  is  controlled  by  a  distribution  of  energy  minima  and  acti¬ 
vation  barriers  associated  with  a  disordered  network.  This  has  been 
revealed,  for  example,  in  the  study  of  oxidation  of  Si  [29,30]  and 
hydrolysis  of  silica  [31  ].  The  present  analyses  of  Li  insertion  in  a-Si 
and,  in  particular,  the  calculated  ring  statistics  suggest  a  simplified 
view  of  the  energy  landscape  governing  the  Li  diffusion  and  mixing 
in  a-Si.  Namely,  the  interstitial  holes  represent  the  deep  sites  (i.e., 
valleys)  in  the  energy  landscape,  and  the  network  incorporation 
sites  correspond  to  the  shallow  sites.  The  deep  and  shallow  sites 
can  be,  respectively,  related  to  the  6-membered  and  3-membered 
rings  that  peak  in  the  ring  distribution  function  at  different  x  val¬ 
ues  (Fig.  3).  While  the  deep  sites  may  play  an  important  role  in 
the  volume  and  voltage  changes  during  lithiation,  the  rate  of  Li 
diffusion  is  likely  to  be  governed  by  the  shallow  sites  with  low 
activation  barriers  between  neighboring  shallow  ones  and  by  their 
percolation  in  the  disordered  network.  A  mechanistic  study  of  the 
dynamic  lithiation  process  would  be  highly  desirable  via  a  multi¬ 
scale  model  that  incorporates  the  above  insights,  e.g.,  by  coupling 
the  atomistic  barriers,  network  topology  and  kinetic  Monte-Carlo 
sampling,  as  implemented  in  the  Si  oxidation  simulation  [30].  This 
is  important  for  understanding  the  electrochemical  alloying,  and 
could  have  implications  on  improving  the  energy  density  and  cycle 
reversibility  of  Li-ion  batteries. 
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